SPACINGS AND PAIR CORRELATIONS FOR FINITE BERNOULLI 

CONVOLUTIONS 



ITAI BENJAMINI AND BORIS SOLOMYAK 

Abstract. We consider finite Bernoulli convolutions with a parameter 1/2 < A < 1 
supported on a discrete point set, generically of size 2^. These sequences are uniformly 
distributed with respect to the infinite Bernoulli convolution measure v\, as TV — > oo. 
Numerical evidence suggests that for a generic A, the distribution of spacings between ap- 
propriately rescaled points is Poissonian. We obtain some partial results in this direction; 
for instance, we show that, on average, the pair correlations do not exhibit attraction or 
repulsion in the limit. On the other hand, for certain algebraic A the behavior is totally 
different. 



1. Introduction 

With original motivation coming from physics [lj, many authors investigated spacings 
and other statistical properties of various number-theoretic sequences, see e.g. [9| [T6 ], [T7 | [TH] 
and references therein. Given a sequence {9 n } equidistributed in the unit interval, the 
nearest-neighbor spacing distribution is obtained by ordering the first elements of the 
sequence 9\ t ^ < 6*2, n < • • • < @n,n, and then defining normalized spacings to be 

$W := N(6 n+ltN - 9 n , N ). 

One can also consider next-to-nearest neighbor spacings or more generally, for any fixed 
£> 1, 

4 N J ■= N(9 n+e>N - 9 n , N ). 

A function Pe(s) is called the limiting distribution function of as N — ► oo if, for 

any interval [a, b] C [0, 1), 

(1.1) lim l # {n<7V: ^ G [o,6]}= / 'p e (s)ds. 
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If the points 6 n are independent uniformly distributed on [0, 1] random variables, we get 
the Poisson model, with 

-(<-!) 

P ^ = ji^Ty 6 ' 3 ' forall ^ L 

Other widely used statistical measures are correlation functions. The pair correlation 
function for the sequence {6 n } is defined by 

R2{s, {9 n },N) = — #j(n,ra) : n^m,n,m<N, \6 njN - miN \ < s|. 

In the Poisson model we have i?2(s, {On}, N) — ► 2s, as N — > oo. 

In this paper, we initiate a study of a similar problem for sequences arising from over- 
lapping iterated function systems and non-standard digit expansions. For A G let 

N-l 

A N (X) = {(1 - A) a n X n : a n G {0, 1}}. 

n=0 

We are interested in the randomness properties of these sets as N — * oo. Note that a 
point in An (A) may have more than one representation; in other words, we may have 
"coincidences." This will happen whenever A is a zero of a polynomial with coefficients 
{—1, 0, 1}, which we call {0, ±1} polynomials for short. In this case we take the point "with 
multiplicity," so that our sequence always has size 2^. Consider the uniform measure 
on An(A); this is a discrete measure, with a point in An(A) having mass 2~ N times the 
number of representations. It is well-known that for all A, the measures z/jj converge 
vaguely to the infinite Bernoulli convolution measure v\, defined as the distribution of the 
random series X^nLo — ^)^ n > w here uj n are i.i.d. Bernoulli random variables taking 
the values 0, 1 with equal probability. Often ±1 are used as digits instead of 0, 1; this is 
just a linear change of variable. We also multiplied by (1 — A), so that the measures are all 
supported on [0,1]. Note that (^,1) is the interesting parameter range for this problem: 
for A < 1/2 the measure v\ is supported on a Cantor set, and An(1/2) is just the set of 
binary rationals in [0,1) with denominators less than 2^; also, for A < 1/2 the ordering 
of A/v(A) agrees with the lexicographical ordering on {0, 1}^. 

The measures u\, for A G L 1), have been much studied, but are still rather mysterious, 
with many unsolved problems, see [12] and references therein. Besides mathematics, they 
come up in physics, information theory, and economics, see e.g. [21 [10]. Briefly, it is 
known that v\ is continuous for all A, it is either purely absolutely continuous or purely 
singular, it is absolutely continuous for a.e. A G (g, 1), and it is singular for Pisot reciprocals 
(see the definition below). It is an open problem whether these are the only exceptions. 
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The only explicitly known A for which v\ is absolutely continuous are Garsia reciprocals 
(see the definition below), which include 2~ l l k , for k > 2. Moreover, the roots of 1/2 are 
the only parameters for which the density is explicitly computable. (When we talk about 
"density of v\ in this paper, we mean the Radon-Nikodym derivative ^p-.) 
Let {£j,7v}f=i be the ordering of Aj^(X), counting with multiplicity, so that 

= £i,jv < &,N < ■■■ < t.2 N ,N = 1 ~ ^ N ■ 

Spacing statistics provide a measure of randomness for the sequences {£n,iv} approximating 
v\. Numerical evidence seems to indicate that, for a typical A, there exist level spacing 
distributions. However, we first need to rescale the sequence. Let 

Cn.JV = Fx(£ n>N ), n = l,...,2 N , where F x (f) = u\(-oo,£). 

These sequences are uniformly distributed modulo one, in the sense that for any interval 
[a, b] C [0, 1], the number of points £ n ,Af m [a, b] divided by 2 N tends to b — a. (This follows 
from the fact that v\ is continuous, hence the cumulative distribution function (CDF) of 
converges to the CDF of v\ everywhere.) 

Conjecture 1. For almost every A € (|, 1) the rescaled sequences ^ Hj n are distributed 
according to the Poisson model. 

We have numerical evidence in support of this conjecture, however, proving it seems 
far beyond our reach. A key difficulty is that very few measures v\ are known explicitly. 
Therefore, we mostly focus on the pair correlation function for the non-rescaled sequences 

R 2 {s,X,2 N ) :=— #|(x,y) : x ^ y, x,y e A N (X), \x - y\ < ^} 

which provide information about the variance of average spacings. A more modest goal 
than proving Conjecture 1 would be to verify the following 

Conjecture 2. For almost every A G (-, 1), there exist c, C > such that 

cs < R 2 (s, A, 2 N ) < Cs for all N and s > 0. 

This would mean that our sequences exhibit no attraction and no repulsion. Our 

main results (Theorems 2.1 and 2.3) makes a step in this direction, but we have to do 
averaging over the parameter A. 
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Note that for certain A the behavior of the spacings is definitely not Poissonian. For 
example, let A = 1/9, where 9 is a Pisot number, that is, an algebraic integer > 1 whose 
conjugates are all less than one in modulus. Such with A G (^, 1), are the only singular 
infinite Bernoulli convolutions known [3J. In this case, Garsia's Lemma [6] says that the 
spacings between distinct points in An(\) are all bounded below by C ■ X N , which implies 
that there are "massive" coincidences, since 2 N — C ■ X N spacings are zeros. This is an 
extreme case of attraction. Something similar happens for all zeros of {0, ±1} polynomials. 

On the other hand, consider the case of "Garsia numbers." We say that 9 is Garsia if 
it is a zero of a monic polynomial with integer coefficients and constant term ±2, whose 
conjugates are all greater than one in modulus. Garsia reciprocals include 2~ 1 / fc for k > 2, 
but also many other numbers (e.g. the reciprocal of the positive root of x 3 — 2x — 2 = 0, 
which is about .5652, or the reciprocal of the positive root of x 3 — x 2 — 2 = 0, which is 
about .5898). Garsia reciprocals A are the only specific (as opposed to generic) A £ (^, 1) 
for which it is known that v\ is absolutely continuous, and moreover, they have bounded 
density. Again by Garsia's Lemma [6], it is known that the spacings between points in 
j4tv(A) are bounded below by C ■ 2~ N , and there are no coincidences. This is an extreme 
case of repulsion. 

It is interesting to compare our situation with that for sequences an d mod 1, for ir- 
rational a and d > 2. It is conjectured that for a badly approximable by rationals the 
spacing distribution is Poissonian, and partial results in this direction are obtained in 
[16\ I18j. in particular, for a.e. a, the pair correlations converge to 2s as N — » oo. On the 
other hand, |16|, |18j show that for some irrational a the behavior of spacings is different. 
In our problem, we prove that for A which are extremely well-approximable by zeros of 
{0,±1} polynomials, the pair correlation function i?2(s,A, 2 N ) exhibits some attraction 
(Theorem 2.2). 

There are several difficulties in pushing our results further; one of them is that the 
distribution of zeros of {0, ±1} polynomials is not well-understood. In contrast with 
[18j . we cannot say anything about higher-order correlations; this is related to hard open 
problems on Bernoulli convolutions. A technical assumption which we have to make in 
some results is that the parameter A belongs to the "transversality interval," that is, an 
interval where {0, ±1} power series cannot have double zeros. Transversality [1511201 fT3], 114] 
has been crucially important in the work on v\. 



SPACINGS FOR FINITE BERNOULLI CONVOLUTIONS 



5 



We note that there are some features of our sequences which are "non-random." In 
particular, some spacings tend to have very high multiplicity, in other words, the distri- 
bution of non-rescaled spacings sometimes exhibit "spikes." This is explained by the fact 
that every spacing g is a zero of a {0,±1} polynomial. For every such polynomial with k 
zero coefficients, the same difference between points in A]y(X) will be achieved 2 k times. 
If g is a very small number, it is likely to be a spacing, and if N — k is small, we will see 
a "spike." This irregularity should dissappear after rescaling by the CDF of v\. 

Further, it turns out that the smallest spacing between points of Ajv(A), for a.e. A in 
a transversality interval, is bounded below by 3~ N N~ x ~ e , for infinitely many N, for any 
e > (follows from Theorem 2.4), which is much larger than (2~ N ) 2 which we would get 
if the points were sampled randomly and independently from the uniform distribution. 

In Figures 1-3 we show the histogram of the distribution of nearest-neighbor spacings 
and also spacings with £ = 2, 3 for the rescaled points {£ n ,7v}; where A = 0.70880447 
was obtained from a random number generator, uniformly in the interval [0.69,0.71]. It 
was chosen so as to be reasonably close to 2 -1 / 2 = 0.7017. . . for which the CDF of the 
infinite Bernoulli convolution measure is known precisely, and this CDF was used for the 
rescaling. The figures show reasonable agreement with the Poisson distribution (indicated 
by the curves). Note that the units on the horizontal axis are 2~ N . More details about 
the computations, Mathematica program, and additional simulations are collected in the 
Appendix at the end of the paper. 
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Figure 1. Nearest-neighbor spacings for rescaled points, A 
.70880447, N = 22. 
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Figure 3. Spacings with I = 3 for rescaled points, A = .70880447, N = 22. 



2. Results 

Theorem 2.1. Let I = [Ao,Ai] C (^,0.668). Then there exists C\ (depending on I) such 
that 

(2.1) V N > 1, Vs>0, i? 2 (s, A, 2*) dX < C lS , 
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Remarks. 1. The appearance of the number 0.668 is due to "transversality;" it is 
explained in the next section. We do not know if the result holds all the way to 1, 
although numerical evidence indicates that things actually get better as A increases. 

2. By Fatou's Lemma, we get J 7 liminfjv_ >00 i?2(s, A, 2 )dX < C\s for s > 0. This 
shows that, on average, our spacing distributions exhibit no attraction, at least in the 
liminf sense. We note that in other problems, for instance, concerning the distribution of 
{na} for irrational a, it is also common to average over the parameter, see [3 Th.3.2]. 

3. An application of Borel-Cantelli shows, for instance, that for any e > 0, for a.e. 
A G (i, 0.668), 

(2.2) Vs > 0, R 2 (s, A, 2 N ) < sN 1+£ for N sufficiently large. 

4. An estimate analogous to (|2.ip follows for the sets An(X) n J, for an arbitrary fixed 
interval J C [0, 1]. Let #{J,N) = #(A N (X) n J) and set 

(2.3) R 2 (s,X,J,N) := ^±— #{( x ,y) : x ± y, x, y € A N {\) D J, \x-y\< ^^y}- 
Then, obviously, 

/ s\J\ \ 2 N 

Note that lim^oo #{J, N)/2 N = u x (J) £ (0,1), so we get an analog of (gZEJ for 
i?2(s, A, J, N), with a constant Ci^ ~ C\\ J\/v\{J) 2 for large N. 

5. If z^a has a density bounded below by rj > on an interval J = [a, b] C [0, 1] , then 

R 2 {s, A, J, AT) < R 2 (s/v, A, J, A 7 '), 

where R 2 (s, A, J,N) is the correlation sum analogous to (|2.3|) for the rescaled points £ n ,N 
in J = [F\(a) , F\(b)] (recall that F\ is the CDF of v\). It is known that the density u\ is 
bounded away from zero on [e, 1 — e] for e > 0, for a.e. A > 2 -1 / 2 , and numerical evidence 
suggests that this is true for a typical A > (y/E - l)/2 « .618 as well. 

The following theorem shows that (|2.2p fails for many points. 

Theorem 2.2. (i) For any Xq which is a zero of a {0, ±1} polynomial, there exists p = 
p(Xo) < 2 such that 

(2.4) |Aiv(Ao)| < C 2P N , 

where \An(Xo)\ is the size of the set without multiplicities. Thus, 

(2.5) J R 2 (0,A ,2 JV )>C7 2 - 1 (2/p) JV -l. 
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(ii) For any interval I C (|, 1) and e > 0, there is an uncountable set £j C /, such that 
for X G Si 

(2.6) Vs > 0, R 2 (s, X,2 N )> 2 Nl ~ s for infinitely many N. 

Remark. One can ask how large is the uncountable set <?/ in terms of Hausdorff 
measure and dimension. Our proof certainly produces a set of Hausdorff dimension zero, 
but it shouldn't be hard to find a gauge function for which the Hausdorff measure is 
positive. We leave this to the interested reader. 

Next we turn to the estimates from below. 

Theorem 2.3. For any Ao G (§,1), e G (0, 1 — Ao), and L > 0, there exists C3 = 
Cs(Ao,e) > such that 

/•Ao+e 

(2.7) V N > 1, V s G (0,L), / R 2 (s,\,2 N )d\ > C 3 s. 

J Ao-e 

Remarks. 1. An analog of (12. 7p holds for R2(s, A, J, N), for any interval J C [0,1], 
with a constant C^j. 

2. If v\ has a density bounded above by M, then the pair correlation for the rescaled 
points satisfies 

R 2 {s,X,2 N ) > R 2 (s/M,X,2 N ). 
It is known [20] that v\ has bounded density for a.e. A G (2 -1 / 2 , 1). 

It is often more convenient (especially for the proofs) to work with 

N-l 

(2.8) A' N (X) = {Y1 a - A ™ : a n 6 {0, 1}} = (1 - A)- X A W (A) 

n=0 

and the corresponding correlation expression 

R' 2 {s,X,2 N ) = ^#{(x,y) : x^y, x,y e A' N (X), \x - y\ < 

Since (1 — A) -1 is uniformly bounded from zero and from infinity in all our theorems, 
proving results for R^is, A, 2^) will immediately imply the desired statements. 

Next we consider the smallest and the largest gaps (spacings) for A' N (X), which we 
denote <?at(A) and Gjy(X) respectively. It is obvious that 5w(A) < (1 — A) _1 2 _iV . For 
Garsia reciprocals A we have <7at(A) > c- 2~ N . It seems likely <7jv(A) = o(2~ N ) for a typical 
A, but we do not have a proof of this. The difficulty, again, is that we do not know much 
about the distribution of the zeros of {0, ±l}-polynomials. 
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We have a lower bound for gjv(A) for a typical A in the "transversality interval." 

Theorem 2.4. Let > be such that *}2 N 3 n on < oo. Then for Lebesgue a.e. A G I = 
[^,0.668] we have 

<7jv(A) > «at /or infinitely many N. 

Concerning the largest gap, we do not have new results; we just note that questions 
related to Gat (A) have been studied by Erdos, Joo and Komornik (though they were 
couched in different terms). It is easy to see that Gat(A) < A^ -1 [U Th.4(a)] (this 
follows by the "greedy algorithm"), and hence Gtv(A) = A^ -1 , since the ordered sequence 
A' N (X) starts with OjA^ -1 . O ne can ask, however, whether this largest gap is realized 
away from the endpoints, as N — > oo. The answer is "yes" for all A G (^,A 9 ), where 
\ g = (V5 -l)/2« 0.618 ... is the golden ratio [U Th.4(b)], in which case, for odd N , 

1 + A 2 + A 4 H h A^- 3 , 1 + A 2 + A 4 + • • • + A^- 3 + A^ 1 

are two consecutive points in A' N {\). Erdos, Joo and Komornik [SJ Th.5] also showed that 
Gat(A) is achieved at most K times, independent of N, and not at all in (e, (1 — A) -1 — e) 
for any e > 0, for sufficiently large N, for all A G (2 -1 / 2 , 1) which are not zeros of {0, ±1} 
polynomials. It is not known what happens between A 9 and 2 -1 / 2 , though it is reasonable 
to expect similar behavior (gaps of o(X N ), except near the endpoints), at least for typical 
A. 

3. Proofs 

We need a result about {0, ±1} polynomials. Actually, it is convenient to state it for 
power series; of course this includes polynomials as a special case. Denote 

oo 

B= {l + ^c n zr n , Cn G {-1,0,1}}. 

n=l 

Definition 3.1. An interval I C (0, 1) is called an interval of transversality for B if for 

every f G B, we have 

(3.1) x€l, f(x) = => f'(x) + 0. 

Proposition 3.2. [19] The interval [i, 0.668] is an interval of transversality for B . 

Transversality is used in the following lemma. The Lebesgue measure on the line will 
be denoted by L. 
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Lemma 3.3. ( \14:\ Lem. 4.2], see also |15j ) Let I be a closed interval of transversality for 
B and g E B. Then there exists C > 0, depending only on I, such that 

(3.2) Vp > 0, £{X € / : \g{\)\ < p} < Cp. 

Recall that we will be working with A' N (X) and H^(s, A, 2 ), see (|2.8p . 

Proof of Theorem \2.1\ Our proof mimics the argument in [13] , but we have to work 
with finite sums instead of power series. Let Oat = {0, 1}^ and tv^ (u>) = Yln=o 0J n^ n for 
oj = ujq . . . un-i € 17 j\f . Note that A' N (X) = it^(Qn)- We will assume that ir\ is 1-to-l 
on rijv, m other words, A is not a zero of a {0, ±1} polynomial, which only excludes a 
countable set from /. Further, let pjy be the uniform (counting) measure on S7jv> with 
each point having equal weight 2~ N . We have, by definition, 

R' 2 (s,X,N)=2 N (p N x p N ){(cu,T)en 2 N : u^t, \ir% (u) - Trf (r)| < s2~ N }. 

By Fubini's Theorem, 
(3.3) 

f R' 2 (s, A, N)dX = 2 N [[ £{XeI : - vrf (r)| < 52^} d^w) ^tv(t), 

where Diag = {(w,w) : w G O}. We have 

N-l 

(3.4) Trf (a,) - vrf (r) = - Tn)A n = ±A fc 5 (A), 

ra=0 

where k = \uAt\ = min{n : uj n ^ r n } < N—l, and g is a {0, ±1} polynomial with constant 
term 1, hence g £ B. By (J33D, |?rf (w) - vrf (r)| < s2 -JV implies |#(A)| < Aq k s2~ N =: p 
for A e I = [A , A x ]. Applying fl32J) yields 

£{A G / : |yrf (w) - vrf (r)| < 82"^} < CX^ k s2- N , 

and we obtain from ()3. 3|) : 

[ R' 2 (s,X,N)dX < Cs ff A |wAr| dp N (uj)dp N (T) 

J I JJnl\T)iag, 



N-l 



= Cs Aq fc (MiV x /iftf){(w, t) : | a; A r| = fc} 

fc=0 

AT-1 „. 

- ^2^ A o 2 < 2Ao^T' 

fc=0 u 

proving (|2.ip . □ 
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Proof of Theorem \2.2l (i) Let Ao be a zero of a {0, ±1} polynomial, that is, 

k 

l + J^c n A n = 0, with c n G {0, ±1}. 

n=l 

We can write c n — u n — v n for some u n , v n G {0, 1}, n — 1, . . . , k. It follows that 1 + 
Y2n=l u n^o = En=i v n^o ' which means that every occurrence (not just at the beginning) 
of the "block" (lu\ . . . u^) in a sequence a$. . . a^-i representing a point in A' n (\q), can be 
replaced by (Oui . . . Uk), without changing the value of the sum, in other words, yielding 
the same point in ^4^y(Ao). It follows that |^4jy(Ao) | is bounded above by the number of 
"words" of length N in the shift of finite type consisting of all 0, 1 sequences which do not 
contain the "forbidden" block (lux . . . u/,), see [5J Ch.2]. Further, it is known (see [8]) that 
the number of words of length in a shift of finite type is not greater than const • p N , 
where p is the topological entropy. The entropy of the full shift on two symbols is 2, and 
every "forbidden word" results in a drop of entropy, see [H Th.4.4.7]. Thus, p < 2, and 
(|2,4p is proved. Now (|2,5p follows from the Cauchy-Schwartz inequality: the 2^ points of 
An(X) (counted with multiplicity) are grouped into at most const ■ p N distinct values, and 
the sum of the squares of multiplicities is equal to 2^(^(0, A , 2^) + 1) (the comes 
from the diagonal). 

(ii) Let Z denote the set of zeros of {0, ±1} polynomials. Observe that Z is dense in 
(|, 1). Indeed, given A G (^, 1), for any k, we can find c\ = 1, cz, ■ ■ ■ , Cf. G {0, 1} such that 
1 - En=i c n A ™ G [°> Xk ) ( e -g- b y the "greedy" algorithm). Then p(x) = 1 - Yl=\ c n^ n 
has derivative p'(x) < —1 and hence p has a zero in [A, A + A fe ). 

We construct points A G /, extremely well approximable by points in Z, to ensure 
that (|2.6p holds. We do this inductively. First find Xq £ Z f] I and find Nq so large that 
#2(0, A , 2^°) > 2^0"" , which is possible by ([23]). Then we find a closed interval I\ around 
Ao, so small that 

#2(1^,2^°) > 2^0"" for A G h, 

which is possible by continuity. We continue this process and find closed intervals I D 
I\ D . . . D Jfc D . . . and positive integers A^o < Aq < . . . < < . . ., so that 

R 2 (2' k , A, 2^*) > 2 N l~* for A G I k , k > 1. 

We are using density of Z and (|2.4p at each step (however, the exponential rate depends 
on the point in Z, so N^+i — Nf, may be quite large). Clearly A G HfcLi wm satisfy 
(|2.6p . Modifying the construction slightly, namely, choosing two distinct points in Z in 
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an interval and two disjoint intervals around them at each stage, we obtain a Cantor set 
of points in / for which (|2.6[) holds. □ 



Proof of Theorem \2.'J[ The claim (|2.7|) will follow if we prove that 

(3.5) Z(s)--= [ ° + R' 2 (s,\,2 N )d\>C f 3 s. 

For any k > 1 we can find c\ = 1, C2, . . . , Cf. G {0, 1} by the greedy algorithm, so that 

k 

l-J> n A£G[0,A fc ). 

n=l 

Choose k so that Aq < e/10. We have k > 4, since Ao > ^. The idea is that for some 
integer £, for N sufficiently large, the pair of points x(X),y(X) G ^4^(A), where 

N-l k N-l 

x(A) = l+ Yl a " A "' = ^c„A-+ £ b n X n , 

n= k+£+l n=l n=fc+£+l 

will satisfy \x{\) - y(X)\ < s/2 N for all A in an interval of size > const • s/2 N contained 
in (A — e, Ao + e), for all choices of a n , b n G {0, 1}. This will contribute > const • (s/2 N ) ■ 
(^N-k—i—i j2 N ) = const • s to the integral in (|3,5p . proving the estimate. Now for the 
details. Consider 

jV-1 

(3.6) /(A) := x(X) - y(X) = 1 - Cl X - c 2 A 2 - . . . - c k X k + ]T (a n - b n )X n . 

n= k+i+l 

We have 

\k+e+i x e+i 

l/(Ao)l<AS+ ^_ < ^ + ^ { A_. 

Assume that £ satisfies 

(1 - Ao)-^/ 1 < 1 
(there will be another condition on £ below). Then 

(3-7) |/(Ao)|<|. 

Now we need to estimate /'(A) for A G (Ao — e, Ao + e). We have 

oo 

(3.8) imi^rf-^a-A)- 2 . 

n=l 

On the other hand, since c\ = 1 and c n > 0, we obtain from (|3.6() : 

/'(A)<-1+ jr nX^<-l+ k -±^X^. 

n =k+£+l \ ' 
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Let Ai = An + e and assume that I G N is such that 



k + e + 2 xk+i i 

(l-Ax)2 Al S 2' 



then 



(3.9) /'(A) < -1/2 for all A £ (An — e, A + e). 

Prom ([321) and §2M it follows that there exists 7 G [An — ^, A + ^] such that 7(7) = 0. 
Then ([31]) implies that 

|/(A)| < S 2- w for A G (7 - (1 - Ai) 2 «2-* 7 + (1 - Ax) 2 ^), 

and we have s2~ N (l — A1) 2 < ^ for large N. Thus, any choice of a n ,b n G {0,1}, with 
n G [A; + ^ + 1, N - 1] D N, contributes at least s2~ N (I - A x ) 2 • 2~ N to the integral J(s), 
and taking all the 4A— fc-«-l cllo i ces in t account yields (|2.7[) . □ 

Proof of Theorem \2.4\ It is clear that 5w(A) = S^o'Con - b n )\ n for some a n , 6 n G {0, 1}, 
hence 5at(A) = A fe |p7v_fc(A)| for some polynomial pat-A; G 23 of degree iV — A; — 1, where 
fc G {0, 1, . . . , N - 2}. If g N (X) < a N , then 

0< \p N -k(X)\ <a N \~ k <a N 2 k . 

By Lemma 13.31 

£{A G / : \p N - k {\)\ < a N 2 k } < Ca N 2 k . 
There are ^^k^ 1 polynomials of degree N — k — 1 in B. Let 

E N = I A G I : 3at(A) < 

It follows that 

AT-1 

£(#iv) < Ca^v ^ 2 ^ • 3 N - k - x < Ca N 3 N . 

k=0 

Recalling that on^ N < 00 we conclude by Borel-Cantelli that a.e 
belong to En for infinitely many N, proving the claim. 



. A G I does not 
□ 
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4. Appendix: numerical experiments 

In Figure 4 we show the histograms for the non-rescaled and rescaled data sets for 
A = 0.6429, which was chosen randomly. For all rescalings we used the density for the 
parameter 2 _1//2 given by 

ax 2 /2, x E [0, b] 

F(x) = < ab(x-b), a; €[6, 1-6] , where a = 1.5^2 + 2, b = V2 - 1. 

1 -a(l -x) 2 /2, x e [1-6,1] 

This "improved" the distribution of spacings even for A not very close to 2 -1 / 2 . In the 
histogram, we used the 50 "bins" from 0.1 • 2~ N to 5-2^. 





rescaled non-rescaled 
Figure 4. Spacings for a "random" A = .6429, N = 22. 

Similar results (though, admittedly, some showing less accurate fit) were obtained for 
more than 20 values of A in the range [0.62, 0.78]. For some values we went up to N = 25, 
but this significantly increased the computing time without noticeable advantage. 

We also tried A close to a Pisot reciprocal and to Garsia reciprocal. In Figure 5 we took 
A = 0.6518 which is close to the Pisot reciprocal 0.651822...; we still see "exponential- 
like" distribution, but there is definite "attraction" visible. For A = 0.652 this effect was 
no longer felt. In Figure 6 we took A = 0.70710678 which is within about 10" 9 to 2" 1 / 2 , 
a Garsia reciprocal. It shows a very dramatic "repulsion." This effect was not felt for 
A = 0.7071. 
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rescaled 



non-rescaled 



Figure 5. Near Pisot: A = .6518, N = 22. 




rescaled 



non-rescaled 



Figure 6. Near Garsia: A = .70710678, N = 22. 



To conclude, we provide the Mathematica 6.0 program used to create these figures. 
We use A =lam= 0.6429, N = 22; I =ell= 1 means that we consider nearest-neighbor 
spacings. We then just change these parameters. Spac and SpacR produce the histograms 
of spacings for non-rescaled and rescaled points, respectively. 
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Needs ["Histograms ' "] 

a = 1.5*Sqrt[2.] +2.; b = Sqrt[2.] - 1. 

F[x_] := Which [x < b, a*x~2/2., x > 1 - b, 

1 - a (1 - x)~2/2., True, a*b~2/2 + a*b*(x - b)] 

lam = 0.6429; N= 22; ell=l; Clear [g] ; g[A_] := Union [lam A, lam A + 1] 

Tru = Nest[g, {0}, N] (1 - lam); 

Vru = Drop [RotateLeft [Pru, ell] - Pru, -ell] 2~N; 

Spac = Histogram [Vru, HistogramRange -> {0, 5*ell}, 

HistogramCategories -> 50, Approximatelntervals -> False] 

Pru = Table [F [Tru [[i]]] , {i, 1, 2~N>] ; 

Vuru = Drop [RotateLeft [Pru, ell] - Pru, -1] 2~N; 

SpacR = Histogram [Vuru, HistogramRange -> {0, 5*ell}, 

HistogramCategories -> 50, Approximatelntervals -> False] 

We compared the histograms with appropriately rescaled Poisson distribution as follows: 

Pois[s_,ell_,N_] : = .l*ell*2~N*s~(ell-l) Exp [-s] / (ell-1) ! 
Plot [Pois[s, ell, n] , {s, 0, 5*ell>, PlotRange -> All] 

Acknowledgment. We are grateful to Rick Kenyon for his ideas at the initial stage 
of this work. Thanks also to Yury Makarychev for some help with simulations. 
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